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An adaptive decision-directed sell erne using- 
circulating filters is described for the estimation of the 
impulse response of fading dispersive channels. Theoretical 
and computer simulation results show that more than 12 dB of 
SITE improvement can be realized. 

Incorporating the above scheme and taking advantage of 
the skewed na ture of the delay power spectrum of troposcatter 
channels s a decision feedback equalization scheme which can 
implement matched filtering with a small loss of the received 
signal energy is considered and its performance evaluated hy 
computer simulation. Also studied is the performance of the 
circulating filters in Viterbi algorithm -Ml sequential decoder. 
Simulation results indicate that a low probability of error 
performance can be obtained. 



CHAPTER 1 


INTRODUCTION 


This thesis considers the problem of adaptive estimation 
and equalization of the fading dispersive channels in the 
contest of high speed binary data transmission. The objective 
of the study is to develop an efficient channel estimation 
technique and to observe how it performs in some representa- 
tive examples of channel equalization schemes. 

1 . 1 Backgr ound 

Most of the communication channels are band limited and 
data transmission over them results in symbol dispersion 
and hence in intersymbol interference. Also corrupting the 
received signals are the thermal and other noises. As a 
result, the received signals bear almost no resemblance to 
the transmitted signals. Efforts are made to minimize these 
effects at the transmitter end by optimizing the pulse shape 
[Gha~50, TS-64] so that it suffers least energy loss in 
passing through the channel (signal shaping) and at the 
receiver end by compensating for the channel dispersion 
(channel equalization) [LSW - 68]. A number of theoretical 
schemes have been proposed (Tit -6 7, AP ~ 70] to permit high 
performance data communication over these channels. A major 
problem faced in practical application of these schemes is 
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that channel characteristics are rarely known in sufficiently 
accurate detail to allow a design that will result in near 
optimum performance. The direct approach of making channel 
measurements and then designing the appropriate receiver is 
impractical in most cases. The technique then adopted is to 
use adaptive receivers that are capable of adjusting auto- 
matically to accommodate changing channel conditions. 

The problem is more acute in tropo-applications where 
the channel is fading dispersive [Bel-63, SBS-66], The 
channel characteristics are in fact the major element to be 
evaluated in the study of any communication scheme for 
troposcatter channels. The performance of any data communi- 
cation system for tropo applications depends to a large 
extent on how well it is able to adapt to channel variations. 

The need for accurate channel characterization in a high 
performance receiver scheme necessiates inclusion of a well 
designed channel estimator in the receiver structure. The 
test-signal and decision-directed methods [1SW-68J are the 
two techniques normally employed for channel estimation. In 
the first case a separate test signal, say a known data 
pattern, sent over the channel is used in extracting infor- 
mation about the channel response. The test signal techniques 
are inaccurate and costly as compared to the decision 
directed techniques. In the decision directed method 
receiver-decisions as to the past transmitted symbols are 
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■used in extracting channel information from the received 

signal and separate test signals are not needed. The 

scheme works well so long as the decision errors occur only 

once in a while. In case of frequent decision errors the 

channel estimates will he inaccurate and equalizer adaption 

will he affected. This may lead to more decision errors. 

In extreme situations, lor,g hursts of errors may result in 

complete loss of channel acquisition, hut in practical 

receivers where probability of error is of the order of 
-5 -6 

10 or 10 this event is really rare. 

The equalization techniques normally employed are linear 
processing and nonlinear or recursive processing of the 
received signal. In linear processing a matched filter 
followed by a transversal filter is used [lSW-68 ].The adjusta- 
ble parameters are the tap gain factors. A number of strate- 
gies like zero forcing method [Luc-65, 1S''/-68J and steepest 
descent algorithm [1SW-68] have been proposed to adaptively 
adjust these tap weights. The performance of these linear 
techniques is significantly inferior to the recursive equali- 
zation techniques like decision feedback equalization [Aus-67, 
GrBS-71] and Yiterbi decoding [Vit-67, Yit~71J. Of 
these, the Yiterbi decoding has been shorn to be optimum 
structure [For-72, hob-71]. 
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The decision feedback equalization is similar to the 
transversal equalizer in that both have a matched filter 
followed by a tapped delay line. However it makes use of 
the previous decisions to coherently subtract the intersymbol 
interference caused by the previously detected symbols. This 
is done by passing the past decisions through the feedback TDL* 
The MF, and forward TDL are used to minimize the effects of 
additive noise and future symbol int erf eronce. 

The Viterbi algorithm is a sequential decoding technique. 
It calculates the likelihood functions of various allowable 
sequences recursively and selects the one with maximum like- 
lihood as the decoded sequence. 

1.2 Pr evi ous Relate d Work 

The 1928 paper by Hyquist [Nyq-28J is a classic paper in 
the data transmission field due to its lasting influence on 
the data system theory and design. The problem he considered 
was to find conditions on the transmitted pulse shape so that 
the samples of the received signal are identical to the 
message sequence in the absence of 1ST and noise. the pulse 
shaping problem is further studied in detail in the subsequent 
papers by Gabon ]_Gab— 46], Chalk |_Cha--50], Tunis and 
Schnidman lTS- 64], DiToro and Steiglitz [DS-65], Ger st and 
Diamond [GD-61] and others and also in books like Bennet and 
Davey [BD-65] and Lucky, Salz and Weldon [LSW-68]. 
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The joint optimization of transmitter and receiver to 
minimize the effects of ISI and noise has been considered 
by Tufts [Tuf-65], Smith [Smi-65J, Berger and Tufts [BT-67], 
Hanster [Han-71], Ericson [Eri-73] and others. 

The work on linear nonre cursive equalizers using 
tapped delay line dates back to Ealmann [Eal-40] who 
proposed use of transversal filters for signal processing. 
Ericson [Erl -71] has shown that for any reasonable criteria a 
matched filter-TDl combination is the espected receiver stru- 
cture. Adaptive versions of linear equalization techniques 
have been given by lucky [Luc-65, Luc-66], Gersho [Ger-69] 
and Proakis and Miller [PM-69J. 

The nonlinear techniques have been studied in two 
different categories. The decision feedback equalization 
has been considered in detail by Austin [Aus-67], Brady 
[Bra-70], George, Bowen and Storey [GBS-71] and Monsen 
[Mon-71, Mon-73, Mon-74 and MR-73]. Brady and Monsen have 
studied the technique for tropo applications. In the second 
categories optimum techniques have been considered, the 
important ones being statistical detection technique des- 
cribed by Abend and Fritchman [AP-70] and Yiterbi decoding 
first described by Yiterbi [Yit-67]. The Viterbi decoding 
has been subsequently studied in detail by Yiterbi [Yit-71], 
Forney [For-72, For-74] and others. A bibliography of work 
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in Viterbi decoding field upto 1973 is contained in 
L'For-73 ] . 

1 • 3 Organization of the The sis 

In tb. is thesis an efficient' channel estimation scheme 
is developed and its performance for two nonlinear equaliza- 
tion schemes,* the decision feedback equalization and Viterbi 
decoding, for fading dispersive channels is considered. In 
Chapter 2 the troposcatter channels are modelled and an 
estimation scheme using circulating filters which is suitable 
for tropo-applications is considered. An expression for noise 
suppression due to circulating filters is derived and simu- 
lation results indicating the estimation performance are 
given. In Chapter 3 a DFE structure for troposcatter channels 
which requires explicit knowledge about the channel response 
is obtained and the performance of the system using circula- 
ting filters for channel estimation is determined by computer 
simulation. In Chapter 4 optimum ML estimator structure 
using Viterbi algorithm is considered. A receiver structure 
using circulating filter estimators is derived and its 
performance is studied by way of computer simulation. 

Chapter 5 concludes the thesis by discussing the results 
obtained and exploring the possibilities of further work 
in the field. In Appendix 1 the details of the simulations 
are given and in Appendix 2 the computer programs for these 
simulations are given. 



CHAPTER 2 


TROPOSCATTER CHANNEL ESTIMATION USING CIRCULATING 

FILTERS 


Signals received through, practical communication channels 
bear little resemblance to the transmitted signals. Efficient 
recovery of transmitted information from such received signals 
depends not only on the choice of signalling waveforms 
transmitted and signal processing techniques employed in the 
receiver but also on the knowledge of the channel chara- 
cteristics. Thus, modelling and estimation of charnel 
characteristics constitute important aspects in the design of 
communi cat ion sy s t ems . 

In this chapter we consider a. scheme of adaptively 
estimating the low pass equivalent complex impulse response 
of a fading time dispersive channel, such as the troposcatter 
channel, over which a binary phase shift keyed (BPSK) signal 
is transmitted. In Section 1 we introduce a fairly general 
model for a fading time dispersive channel . In Section 2 
we formulate the channel estimation problem and indicate that 
the opt. solution implied by such an approach is not 
practical. Assuming the availability of error free decisions 
on past symbols we then illustrate in Section 3 the possi- 
bility of using a circulating filter for estimating the channel 
impulse response. In Section 4, we relax the assumption on 
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intersymbol interference and. discuss a scheme of estimation 
with circulating filters to obtain essentially noise free 
estimates of the channel impulse response. We also obtain 
an expression for noise reduction attainable in the estima- 
tor. In Section 5 we present some simulation results rela- 
ting to convergence of the estimates in the presence of 
symbol decision errors and time-jitter. 

2.1 A Model for Fadin g _D is persiv e Cha nnels 

Following Eailath [Kai-60] and Bello [Bel -63, Bel -69] 
we assume that a fading dispersive communication channel 
can be modelled as a, random linear time— vary ing filter 
with low pass equivalent, possibly complex, impulse response 
h(t,T) which is the response at time t due to an impulse 
applied x sec. earlier. A troposcatter channel which is of 
major interest to this thesis provides a practical example 
of a fading dispersive channel. Empirical evidence as well 
as central limit theorem based arguments [Bel -69, SBS-66] 
suggest that the impulse response h(t,T) of a tropo channel 
may adequately be considered as a zero mean complex 
Graussian process. Further under the usually realistic assum- 
ption of wide sense stationarity and uncorrelated 
scattering such a channel is completely characterized by 
its path-gain correlation function, [Bel-69, SBS-66], 
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E[h(t,T)h* (t+s,r))] = Q(s,t) S(t)-t) 

Moreover, when the fading rate is small compared to the 
information transmission rate the channel may he considered 
to he quasi time-invariant. The delay power spectrum 
Q(0,r), often denoted by Q(r), of the channel is then 

E[h(T)h*(n)] = Q(t ) S(t)-t) 

and represents power spectrum of the channel transfer function 
in frequency domain. 

The delay power spectrum of a typical tropo-link, shown 
in Tig. 2.1, can he reasonably represented by the empirical 
expression 

Q(t) = A't e“ aT t > 0 

2.1.1 BPSIC Transmission over Tropo-channel 

The baseband model of data system is shown in Pig. 2.2. 

The low -pass equivalent impulse response h(t) of the 
tropo-channel is finite, h(t) = 0 for t >_ nT and has n 
chips 

h( t ) = h Q (t) + h-^t-T) + h 2 (t-2T)+ ... h^Ct-n^XT) 

where the chip h^(t) ; i = 0,1, .... (n-1) is nonzero over 
the int erval (0,1], 





Nanoseconds 


. 2.1 Delay power spectra for the short and 
long channels 


Nois e 



,g. 2.2 Baseband model for a communications channel 
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The channel pulse response or composite channel impulse 
response h c (t) is defined as 

h c (t) = h(t) * p(t) 

where p(t) is unit pulse of duration T. 

The input to the channel is a sequence of randomly 
selected equally likely pulses of amplitude + 1. The channel 
output is 

r(t) = £ a, h(t-kT) + n(t) 
k K 

or, r(t) - £ h 1 _h c (t-kT) + n(t) 

k K 

where b^ is input sequence of + 1 impulses and n(t) is zeromean 
complex white Gaussian noise of variance o . 

2.2 Chan nel Estimation Pro blem 

The channel estimation problem is given the received 
or observed signal 

r(t) = £ a v h(t--kT) + n(t ) 
k * 

determine the channel impulse response h(t). 

Optimum solution to the channel estimation problem can 
be obtained by techniques like MAP estimation and Zalman 
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filtering. These methods are considered in detail in 
VanTrees, Part III. Since both the input sequence a k and 
channel response h(t) are unknown, the state equation formula- 
tion of the estimation problem suitable for Kalman filtering 
is nonlinear. The Kalman filter is unable to treat nonlinear 
equations per se and quasi-linearization assumption is to be 
made. The estimator structure is quite complex and the 
technique is not suitable for adaptive channel estimation. 

However, if error free decisions about symbols a k are 
available simpler channel estimation schemes can be derived. 
This we do using circulating filters in the next section. 

2 . 3 Channel Estimation Using Girculating Kilters 

In this section we consider estimation of a non- 
dispersive channel. The received signal in the interval k 
over such a channel is 

r k (t) = a k h(t) + n k (t) 

In the absence of noise 

r k (t) = a k h(t ) 

and h(t) can be determined by correlating r k (t) with a k 

h(t) = a k r k (t ) 
as 




In the general case a noisy estimate h' (t) of li(t) is 
obts.in.ed 

h* (t) = a k r k (t) = a k [a k h(t) + n k (t)j 

= h(t) + e k (t) 

where e k (t) are uncorrelated zero-mean random variables 

2 

with a variance e . 

Essentially noiseless estimates £.(t) ca,n be obtained by 
averaging h'(t) over several, bits. This is done in a 
circulating filter by coherently s.dding its input in succe- 
ssive intervals. The circulating filter scheme is shown in 
Fig. 2.3. For a feedback path gain f, the signal at the 
output of the circulating filter is 

z (t ) = h' ( t ) + fh' (t ) + f %i' (t ) + ... + f^h' (t ) 

= x(t) + y (t ) 

where, for sufficiently large lc, 

x ( t ) = h ( t ) [1+ f + f ^ + ... + f ^ j = 

P ]£ 

and y (t ) = e Q + fe^ + f e^ + .... + f e k . 

p 

The variance cp of y(t) is given by 
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The channel estimate is 

h(t) = (1-f) z(t) 

= li(t) + (1-f) y(t) 

The noise power is suppressed by a factor 

(1-f) 2 <p 2 _ l^f 

e 2 1+f 

iif 

and the corrsponding SER improvement is by a factor j-. 

For a feedback ga,in of 0.9, a substantial noise suppression 
of 12.8 dB is achieved. Thus a channel estimator using 
circulating filters is capable of giving essentially noiseless 
estimates in the absence of ISI. 

In the next section we extend the concepts of channel 
estimation using circulating filters to channels with ISI. 

We consider a more realistic situation -where the input 
symbols are not known but their almost error-free estimates 
ar e avai lab le . 

2.4 Estimation of Dispersive Channels 

In this section we consider estimation of a channel 
whose pulse response h c (t) extends to (n+1) symbols durations. 
The (n-t-l) chips of the pulse response 

= l h.(t - iT) 
i=0 
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are estimated by (n+1) circulating filters one for each 
chip. The schematic diagram of the channel estimation scheme 
is shown in Fig. 2.4. 

th 

The received signal in the k interval is 

n 

■ J 0 Vi h i (t> + 

To A 

In estimating 3 “ chip h. (t)/ 3 = 0 , 1 , .... n in the 

J 

3 filter the interference contribution due to chips 

h?(t), i £ 3 is removed from r k (t) using the symbol estimates 

b-^_ i ; i i 3 and the estimates h?(t) ; i £ 3 corresponding 

to interval (k- 1 ). The resulting signal is correlated with 

the symbol estimate b, • to generate a noisy estimate 

*“* J 

&?' (t) of the chip h?(t). 

J J 

h? (t) = [r^(t) ~ £ h^ ( t ) ]b^_j 3 = 0,1,2, ... n. 

= h?(t) + e^ k (t) 

where the noise torm e--, (t) has components due to additive 
noise n-^t) and the remanent 181 . The e^^.( t ) can be consi- 
dered zero moan uncorrelated random variables with varia- 
2 

ce e r 

The circulating filter then operates upon h° (t) to 
generate the signal z-(t) from which an essentially noiseless 

t] 




Delay 
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estimate ho(t) of h?(t) is obtained. 
2 3 


h-°(t) 


,(t) = + y^) 


where 


y 3 (t) “ °jo + fe jl + f ®3 2 + ** 


with a variance of 9 . 


So that 


k° ( t) 


(1 “ V z 3 Ct) 


= h°(t) + (1-f..) y^(t) 3 = 0,1# •••• n « 

The resulting noise suppression is by a factor 


(1 


_p \ 2 2 

iiLJi 

2 

£ j 


Tlic noise suppressing property of the circulating 
filter in the presence of ISI is considered next. We consider 
a case where the ISI is due to past two input symbols. 

Three circulating filters are used to estimate three segments 
h^(t) § i = 0,1, ' nd of the channel pulse response, 

h c (t). The system can be represented in matrix form as 
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! x ° (k) ; i f ° <\\-i - f ' 2 \\-2 ] 

H 

M 

O 

. 


^(k) 

- | f oVk-l U - f 2 b k-l\-2 

x 1 (k-l) 

( 

! 

Vl 

x 2 (k) 

[ -f o b k b k- 2 " f i b k-l b k -2 f 2 

Xp(k-l) 

t 3 

i 

1 

\ 

_^k-2 _ 


whore r k = + n k 

= feedback gain of the ith circulating filter 
i = 0,1,2 
q = 1 - fi i = 0,1,2 

xCk) = filter output at tine k; i = 0,1,2. 

Or in matrix notation form - 


k(k) - F t(k-l) + I r k . 


In case of no noise 

X(k) = 7 I(k-l) + B s k 
¥e define error matrix as 

S £ (k) i 2(k) - S(k) 
where x_(k) = Xi(k) - x(k) , 

c. j J X 

So that 

X £ (k) = 7 ~ £ (k-l) + n k S 


¥e define error variance matrix as 
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p(k) = E[:: e (k) Z £ (k)] 

^ ,, m 

Substituting for X (k) and X_(k) 

o o 

P(k) = S [ i PS £ (k-l) + n k B J [F - £ (k-l) + n k B] T ] 

As source and noise statistics are indopondent 

P(k) = E[P X £ (k-1) X^(k-1)P T ] + a 2 3 3 T 

2 

whore a = noiso variance. 

= F P ( k“l ) P T + a 2 B 

This after matrix multiplication and noise normalization and 
using the properties that the matrix P(k) is symmetric and 
[t> k ] = 1 for all k can bo written in the form of (1). 

VJe assume that the steady state solution to the above 
system of equations is of the form 



That P. • of above form is really a solution to (1) is seen by 
substituting (2) in (1). Thus P qo under steady state 
assumption is 



2i 




? 

t 

J 

! 

i 

i 

I 

t 

i 
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p ooW = + + l'2 Z ? 22 (k-l) 


2£ o t i b k b k-l ? 10^ k “ 1 ' ) “ 2l o f 2 b k b k-2 :p 20^ k “ 1 ^ 
+ 2f i f 2 b k-l b k-2 P 21^ k-1 ^ + 1 


or, 

o ( l“f -1 ) 

p oo (k) = f 0 ^ + 


o' v o 12 


+f 0 -ry + n-f-jTCf'o+^+fpiT 


d-f 2 ) 2 

+ Ti-i^TTf 0 +f x +f 2 -i ) + 1 

= rr-f 0 )(f 0 H-f 1+ £ 2 -i) t f o + (1-f o )(1 - f i ) 

+ (X-f 0 )(l-f 2 ) + ( 1-f 0 ) (f 0 +f 2 -l ) ] 



p ooM 


Tlie error variance coefficients P_. . are measure of noise 

u 

2 

power ,* since for i = j , qn = P... flic noise suppression 

J j J 


due to j 


-'’ ck circulating - filter is by a factor 


(1 - f,) 2 x “ - 

3 (l-f.)(f +f 1+ f 2 -l) ' 

1-f. 0 

~ f -i-f-i J rf o"l 
o 1 2 


0,1,2. 
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The corresponding SES improvement is 


SITS.. 


mp 


= 10 log. 


f o+f i+f 2 

io 


dB ,* j 


0 , 1 , 2 . 


For f j = 0.9 3 = 0,1,2 / the SITS improvement when there 

is IS I due to past two symbols is 


SSS. = 10 log loL 17] = 12.3 dB. 

A similar result when there is ISI due to past three 
symbols is 


f +• £ -hf -hf —2 

SNR imp = 10 log 10 i-~°~TZ TIT' 3 —] d3 f 2 = 0,1, 2,3. 


Bor fj = 0.9 ,* i = 0,1, 2,3 the SITS improvement is 12.2 dB. 

These results indicate that the penalty is only about 
0.5 dB when ISI is due to past 2 or' 3 symbols and still 
substantial SITS improvement of more than 12 dB is achieved 
by using circulating filters. 


For most practical channels, especially for tropo- 
applications, the estimator must be capable of estimating any 
unknown channel and tracking any variations therein in a 
finite number of steps. That this stability criterion is 
satisfied by circulating filter estimators is demonstrated 
by considering filter equations with any arbitrary initial 
conditions X^(0) 
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X-^l) = F i l( 0) + B r 0 
The system with correct initial conditions is 

1(1) = F X (0) + f r 0 
It is shown that the error matriz 

Z(k) = I (k) - f-j^k) 

tends to 0 as k becomes larger and larger. Detailed deriva- 
tions will be given in some future work. 

In our discussion we have assumed that error-free 
estimates b^. are available. Since the estimation technique 
is decision - directed its performance will be affected by 
frequent decision errors. However, in practical situations 
the error event is very rare and the estimator performance is 
only negligibly, if at all, degraded. Since all practical 
receivers make decisions after a delay of t sec. (decision 
delay), same delay will appear in estimates also. In all 
situations of interest the channel varies slowly as compared 
to data rate and the system is not signif icantly affected by 
using' delayed estimates. 

2*5 Simulati on Results 

The channel estimation scheme described in this chapter 
was simulated on a digital computer (IBM-7044). The channel 
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response estimated had IS I due to past three symbols and was 
corrupted by white G-aussian noise. The feedback: path-gains 
of all circulating filters were 0.9 and four channel samples 
per ba^kl were estimated. The performance of the estimator 
was studied in the following contents 

a) SNR Improvement 

The SNR improvement due to circulating filters was about 
12 dB for all the four circulating filters. This closely 
follows theoretical results. 

b) Convergence 

The estimates were found to converge to correct values 
within 50 data bits when the starting values were zero. It 
took somewhat more tine in attaining the 12 dB SNR improve- 
ment value. For trailing segments (segments with lower SNR's) 
the 12 dB value was reached earlier than for segments with 
higher SNR's. 

c) Decision Errors 

The effect of erroneous decisions on channel estimates 
was studied by deliberately introducing bursts of decision 
errors. It was observed that bursts of upto 5 errors do not 
significantly affect the channel estimates. But with a 
burst of 10 errors the channel acquisition was completely lost. 
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&) Effects of Time-jitter 

‘Time- jit ters of T/4 and T/2 were introduced in the 
received signal and the effect on channel estimates was 
studied. With +T/4 jitter all channel samples are 
correctly estimated except for every fourth sample which 
now depends on data sequence. With -T/4 time jitter, 
first sample of all T sec. segments of channel response 
are erroneous. With time jitters of T/2 sec only two out 
of four samples per T sec. arc- correctly estimated. 



CHAPTER 5 


DFE RECEIVER USING CIRCULATING FILTERS FOR TROPOS CATTER 


CHANNELS 


3.1 Rec ei vers for Tropos catter Channels 

The fading dispersive character of tropos catter channels 
results in the received signal fading and its dispersion over 
several symbol durations. Fading reduces protection from 
thermal noise, whereas dispersion leads to intersymbol 
interference. The channel characteristics are in fact a 
major element to be evaluated in considering any communica- 
tion system for tropos catter channels. In this context 
the channel estimator assumes greater significance in the 
receiver configurations. 

The performance of digital communication systems over 
tropos catter channels is thus degraded predominantly by ISI 
and fading besides additive noise. The degrading effects of 
noise and fading can be reduced by increasing the transmitted 
signal power and by explicit diversity techniques. But this 
way the ISI degradations cannot be minimized since, for 
linear channels, the ISI is a function of channel dispersion 
only and is independent of signal to noise ratio. Special 
techniques, known as channel equalization,, are to be incorpora- 
ted in the receiver to combat ISI. The equalizer is made 
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adaptive to accommodate any changes in the channel 
characteristics. Thus* an adaptive equalizer is a special 
feature of receivers for trope channel applications. 

The channel equalization techniques can he broadly- 
categorized as linear and nonlinear techniques. The linear 
techniques like matched filtering followed by transversal 
filter are straightforward and economic but their performance 
is markedly inferior to nonlinear techniques like decision 
feedback equalization [G-BS-71] and Viterbi algorithm [Vit-71]. 
Especially, the Viterbi algorithm technique has been shown 
to be optimum Ml estimator [For-72, Eob-71]. 

Of the various equalization techniques proposed for 
troposcatter channels, the decision feedback equalization 
and Viterbi algorithm have assumed greater importance; 
the first one because of its practicability and fairly good 
performance and second one because of its optimality. Though 
VA receiver is optimum, there are certain difficulties in its 
hardware implementation like realization of a whitened 
matched filter, increasing complexity for largor dispersions 
of the channels and speed of operation. 

In the present chapter wo study the performance of the 
DFE receivers using circulating filters for tr op o charnel 
applications and the study of the performance of circulating 
filters in a VA receiver for unknown dispersive channels is 



29 


done in the next chapter. 

3.2 DPS Recei ver 

The decision feedback equalization technique was first 
proposed by Austin [Aus-67]. The opt inrun receiver derived 
by hin is very conplex and cannot be inplenented because 
it requires correlation of the output of an infinite length 
transversal filter with all possible data sequences and also 
large number of exponential or s. Practical versions were 
proposed in subsequent papers by Brady [Bra-70], Monsen 
[Mon-71, Mon-73, MR-73, Mon-74], George et al [GBS-71] and 
others. Brady and Monsen have considered DIE receivers for 
troposcatter channels. Brady developed a system which is 
MMSE processor and uses an adaption algorithm based on 
steepest descent method as given by Widrow [Wid-661. The 
receiver structure developed by Monsen is capable of handling 
data at higher rates and adapts to channel variations without 
any explicit knowledge about the channel impulse response. 

The decision feedback equalizer works on the principle 
that the receiver decisions can be used in eliminating past 
symbol interference. The assumption involved is that the 
receiver decisions are- essentially correct. This is a 
valid assumption for low bit error probabilities of the 
order of (10“'’ or 10”^) encountered in practical systems. 
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L receiver scheme using decision feedback equalization 
is shown in Fig. 3.1. The receiver decisions are passed 
through, the feedback filter and its tap gains are adjusted 
to minimize the past digit interference. The renanant ISI 
due to the future symbols is minimized in the forward filter. 
The forward filter is also in the form of a tapped delay 
line and its tap gains are adjusted to minimise future symbol 
interference. 

The received signal r(t) is passed through a filter 
h(-t) matched to the channel li(t) and the output is sampled 
to give the signal 

oo 

yjj. = / r(t) h(t-kT)dt 

— oo 

The forward filter is characterized by tap gains 

i = 0,1, IT and the feedback filter by tap gains 

dj f j =1,2, II* 

IT 

°(t) = l o 6(t-i'T) 
i=0 1 

* M 

d(t) = l a. 

3=1 3 

Then the equalizer output is 

M 

N 

b,' = y o. y . - y d. b, . 

k ±io 1 k+1 )k 1 k - 3 



Tapped delay line 



Tapped delay line 
Decision feedback equalizer 
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The equalizer output is converted, to a decision t>^ 
with, the help of a clipping circuit with zero bias. In 
other words, the equalizer output which ferns sufficient 
statistics for making decision about b^ is checked against 

a zero threshold to give decision b^. 

1 

-1 

Ihe tap gains must be adjusted to minimize the probabi- 
lity that b^. 4 b^.. However, in the absence of analytical 
expressions for error probability in terms of tap gains, 
direct minimization of error probability is not possible. 

Ihe tap gains are then selected to minimize the output mean 

2 

square error, BJ_e^] where 

e k - b i - \ 

The problem of adaptively adjusting the tap gains in 
MMSE sense is considered by George et al [GBS-71] ar\d 
Monsen [Mon-71, Mon-74]. Monsen combined the matched filter 
and the forward TDL into a single forward filter 
The tap weights, c^, are then adjusted to minimize the 
effects of additive noise and future symbol interference. 

The forward filter is tapped at intervals on the order of 
the bandwidth reciprocal because it is approximating a 
matched filter as part of its response. For systems the 
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tap spacing is T/2. The forward filter is then given by 

N 

Y (t ) = y C 6(t - il/2) 
i=0 1 

This T/ 2 tap spacing results in aone perfomance degrada- 
tion due to inherent tine quantization. Errors in tine 
synchronization result in further degradation. Accurate 
synthesis of the forward filter in Monsen's structure is in 
fact not simple. 

Monsen [ Mon-74 jhcsty wqy of simulation shown that the 
DEE receiver for tropochannels performs about 5 dB worse 
than the optimum one shot receiver. The degradation is 
attributed to reasons like only partial ISI cancellation, 
use of a finite number of taps and non-optimum sampling tine 
locations. The last two factors result in inefficient reali- 
zation of the forward filter. It has been pointed out in 
[Mon-74] that the remanent ISI penalty is not more than 1 dB. 
Then the major contribution of a.bout 4 dB to the performance 
loss is from inefficient synthesis of the forward filter. 

In Monsen's scheme [Mon-71, Mon-74] a separate algorithm 
is used for the adaption of the equalizer to the channel 
variations. The adaption is also possible using the 
ch an nel impulse response estimates required for matched 
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filtering of the received signals and a separate of 
adaption algorithm is not needed. 

3.3 DIE Re ceiver using Circulating Filters for Tropo c hannels 

The Monsen's scheme of DEE receiver for troposcatter 
channels [Mon-71 , Mon-74 ] s as mentioned in the proceeding 
section, suffers from a performance loss due to inefficient 
forward filtering. The degradation is substantial and it is 
desirable to use a technique in which the forward filter is 
either better synthesized or it is totally done away by 
using some special property of the channel. In this section 
a receiver structure is considered which does not use a 
forward filter and utilizes channel estimates obtained by 
using circulating filter estimator of chapter 2, for both 
matched filtering and decision feedback equalization. 

In deriving the receiver structure advantage is taken 
of the asymmetric nature of the power impulse response of 
the troposcatter channels. The delay- power spectra of 
typical tropo -links are shown in Fig. 3.2. Fig. 3.3 depicts 
the received signal energy as a function of tine. 

These plots show that the delay power spectrum is highly 
skewed and almost all signal energy is received in the 
first few symbol intervals. As much as 80 percent 

or even more of the signal energy of a pulse transmitted 
over a tropo channel is received in the first T sec. 





Signal 



Figure 3.3; Received. Signal Energy over a Troposcatter Channel 
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Then, if the decision about the transmitted symbol in the 

th th 

k interval is taken at the end of the k interval, there 

will be no contribution due to the symbol transmitted in 

(k+1) interval but sufficient signal energy would be 

available for decisions. Thus by deciding about the transi 

mitted bit at the end of T sec, the symbol duration, the 

future symbol interference is eliminated. The equalizer 

has to overcome the past symbol interference only for whicJq 

the feedback filter alone is sufficient. We do lose some 

signal energy in making decisions at the end of T sec. but 

the resulting degradation is small as compared to that 

incurred with late decisions when a forward filter is a must. 

3.3.1 Receiver Structure 

A schematic representation of the adaptive DFE receiver 
with the above mentioned modifications is depicted in Fig. 

3.4. The receiver has four b_ocks j a channel estimator, a 
decision feedback equalizer, a matched filter and a threshold 
detector. The received signal is delayed by T sec. in the 
channel estimator and then the information about the channel 

response is extracted from it using circulating filters. The 

• • _ 

channel estimates h? are available after a delay of T sec. but for 

J 

normal data rates of information transmission over slowly 
fading tropo channels, this delay is immaterial. The 




/ 
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feedback filter is in the fom of a tapped delay line with 
n taps. The tap gains are set equal to the chip estimates 
of the composite channel impulse response. The receiver 
decisions are passed through the feedback filter. At time 
k the filter output is 

u k< t} = 1 Vj h 3 (t) 

which corresponds to the past symbol interference in the 
received signal 


v k (t) = ^ Vj h j (t) 

The received signal r k is 

* k (t) = l b k _, h?(t) + n v (t) 
j =0 J 0 

= \ b°(t) + v(t) + n k (t) 

Then, the output of the equalizer is 

n 

b^(t) = b k h°(t) + l [b k _. h?(t) - b k _. h?(t)] + n k (t) 

3 

= b k h°(t) + w(t) + n k (t) 
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where the remnant ISI w(t) is 

w(t) = v(t) - u(t) 

It will be quite snail when receiver decisions are correct 
and the channel estimates are fairly accurate. The equalizer 
output is matched filtered by correla.ting it with the locally 

'*Q 

generated chip h . The sampled output g^ of the matched 

filter at the end of kth interval, is sufficient statistics 

for making the decision about the transmitted symbol b^ and 

a decision b^ as to b^ is taken by the hypothesis test. 

1 

^ 00 • * 

7 0 where g^ = / b^t) h c (t). 

—1 — 00 

3.5.2 Diversity Deception 

Diversity reception is also incorporated in the system. 
Maximal ratio combining is used in combining the received 
signals on different branches. Per bth order diversity 
the received signals on 1 branches in che kth interval are 
[r^(t) ,• i = 1,2, ..... L]. Each of the received signals is 
acted upon by seperate channel estimators, equalizers and 
matched filters matched to corresponding channels. The 
sampled HP outputs at the end of kth interval are 
[g^, i = 1,2, ..... 1J. The sufficient statistics as to the 
transmitted symbol is then 
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and decisions 


L 



b^» are nade according to the rule 



0 . 


3.4 P erf oraan ce o f the DI E Receiver for Trop es patter Channels 

The performance of adaptive DFE receiver is limited by 
a self noise due to the adaption process and by intersynbol 
interference. The first effect can be nade reasonably 
snail if the channel is slow fading, i.e. if channel varia- 
tions occur at a rale much less than the data rate. Under 
this condition the behaviour is limited by ISI. I:i addition, 
error-propagation also affects the performance. Due to 
the nonlinear receiver structure, it is not easy to get 
closed form, expressions for probability of error. The 
system performance is then often expressed either in terns 
of performance bounds or simulation results. 

The performance bounds are obtained in terns of the per- 
formance of a one-shot receiver. The one-shot receiver is 
defined as the receiver that detects only a single digit. 

This idealization ignores ISI and provides an error proba- 
bility bound for sequential reception of digits. For no 
diversity case it is simply a matched filter receiver. 1 
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dimensional extension of the W receiver gives one-shot 
receiver for 1th order diversity reception. 

The probability of digit error for one shot receiver 
is given by 

p(y) = i erf c 0/y) (!) 

where y is signal energy to noise power ratio E/N . IT 
is single sided spectral density of the additive Gaussian 
noise process. The received signal energy is given by 

E = J |H C (f )! 2 df 

— CD 

C / \ 

where H (f ) is fourier transform of the composite channel 
impulse response h c (t). 

Since the channel is random, y is a random variable. 

An important performance measure is the probability of error 
averaged over y. For a given average E/iT 0 average probability 
of error is determined by integrating (1) over the probability 
density function for y. This density can be found by consi- 
dering a sampled representation of transfer functions H c (f). 
This is done in troposcatter channel context by Monsen 
[Mo n-73] and an expression for average probability of error 
P(|i), for an independent branch Lth- order diversity system 
with an average E/N = |i(> 1 ) is obtained in the form of a 
convergent series 
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r( i+0.5) 


(l+li) 


i+0.5. 


3.4.1 Error Propagation Effects 

In the receiver configuration described in the 
previous section, the receiver estimates of the transnitted 
bits are used in the removal of the inter symbol interference. 
These decisions are also used in the channel response esti- 
mation. If there is any erroneous decision made, it will be 
propagated back by the feedback filter of the equalizer and 
the intersymbol interference will be only partially cancelled- 
Under such circumstances , there is a higher probability of 
receiver making further errors in its decisions. Error propa- 
gation in channel estimator affects the channel estimates. 
Incorrect channel estimates result in incomplete ISI cance- 
llation and inadequate matched filtering thereby further 
increase the error paprobability. This, any error in 
decision by the receiver increases the cn cnees uf further 

errors in near future and the decision errors tend to occur 
in bursts. However, for DFE receivers, it has been pointed 
out in LGrBS~71] that the receiver is able to recover from 
error conditions and long bursts of errors occur only rarely. 
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The performance of the DEE receiver under error propaga- 
tion conditions has been studied by Austin [Aus-67j» Monsen 
[MR-73 and i v Ion-74] and Gantoni end Butler. [OB-76] using 
two different approaches. In the first approach adopted by 
Austin and Monsen the error propagation is modelled as a 
discrete Markov chain process and bounds on mean recovery 
tine and bit error probability have been obtained. But it has 
been pointed out [CB-76] that the statistics of the noise are 
critical to the development of a Markov model and that the 
method is suitable for error probability calculations for short 
feedback paths only. The approach adopted by Gantoni and 
Butler is based on the theory of success runs. Using this 
technique, they have derived bounds for the probability of 
recovery from error in a given number of steps and for mean 
recovery time which do not depend on the statistics of the 
noise and particular selection ef the channel. 

3.5 Simulation Results 

A digital communication system for dispersive fading 
channels using the DEE receiver described in this chapter 
was simulated on digital computing systems - IBM - 7044 and 
DEO-IO. Details of simulation programs are given in Appendix 
I. Tlie simulation results are given below. 

a) Average Probability of Error - ho Diversity 

The average probability of error characteristic of the 
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system as a function of average signal to noise power ratio 
is depicted in the plot of Fig. 3.5. Also shown is the 
average error rate performance of an optimum receiver for a 
non-dispersive flat fading channel. 

The average prob. of error for the system simulated 
was obtained by averagaug tnc results over one million data, 
bits and 20 different realizations of the channel. These 
results are tabulated in Table 1. 

b) Average Probability of Error -Qua. d Diversity 

A diversity reception of the order of four was considered 
and results for 20 different realizations of the channel over 
one million data, symbols were obtained for four SITE, values. 
Errors were- obtained only for the lowest 31TR va.lue considered; 
i.e. 6 dB when 3 errors occur eel. oor higher Slip’s no 
decision errors were observed. The results core given in 
Table 2. 

c) Error Propagation 

The error propagation effects on the performance of the 
receiver were studied by deliberately introducing bursts of 
'errors in the system. The results are given in Table 3. 
Bursts of upto 5 errors resulted in no extra errors. Bursts 
of 6 error resulted in 2 additional decision errors, burst 




Table 1 Probability of error performance - No diversity 
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Table 3 Error propagation effects 


Burst length. 

No . of error s 

1 

0 

2 

■ 0 

3 

0 

4 

0 

5 

0 

6 

2 

7 

5 

8 

68 

Number of errors 
results are with. 

in 100 bits after burst occurance. The 
quad -diversity and SNR = 3 dB. 

Table 4 

Time-jitter effects 

Time “j it ter 
(.in units of T) 

No . of errors in 

104 bits £ 

0 

0 

+1/4 

0 

-1/ . 

1 

+1/2 

374 

-1/2 

261 


* 


Average SM = 3 dB Quad diversity 
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of 7 error- resulted in 5 additions.! errors end larger 
bursts resulted in complete loss of data acquisition. The 
average SifR at the receiver input was 3 d3. 

d) Tine-jitter 

The receiver performance was studied under the effects of 
tine-jitter by introducing tine-jitters of + T/4 and + T/2 in 
the received signal. The results are given in Table 4. 

The results are for an SM of 3 dB. 



VITERBI ALGORITHM RECEIVER USING CIRCULATING 

FILTERS 

In most of the practical data communi cation systems 
symbol-by-symbol decisions are made mainly duo to the simpli- 
city of the operations involved. However, since better 
performance can be achieved with sequential decoding it i*s 
desirable that the receiver decisions be based on the entire 
received sequence. Unfortunately, the straightforward 
statistically optimum receiver has a complexity exponentially 
increasing with data length [Hel-SO] end is not practical. 

The Viterbi algorithm |_Vit-67j provides a way out and gives 
maximum likelihood estimates of the transmitted sequence 
with greatly reduced efforts |_For-72]. The algorithm is 
based on the fact that practical channels can be represented 
as finite state machines. Since a finite state machine has 
finite number of state— transitions only chese many transitions 
are to be considered by the receiver. Thus the Viterbi 
algorithm is an optimum receiver with fixed complexity which 
bases its decisions on the entire received sequence. 

The algorithm was first proposed by Viterbi in 1967 
[Vit-67] for decoding convolutional codes. Later on the 
method was shown to be a variation of the dynamic programming 



problem by Onura [Omu-69] and to bo a maximum likelihood 
estimator by Forney [For-72] and Kobayashi [Eob-71]. The 
technique has been studied in detail by Viterbi [Vit-71] and 
Forney [For-73] in the context of convolutional codes. In 
another paper j_For-73], Forney presents a general discussion 
on the application of the algorithm to decision problems. In 
more recent papers like those of Hager and Proakis [MP-73J, 
Qureshi and Newhall [QE-73] and Hager and Falconer 
[MP-76A, MF-76BJ attention is directed towards increasing 
the practicability of the 7 A receivers by presenting it with 
a limited memory channel. 

In this chapter a Viterbi Algorithm receiver using 
circulating filters for unknown dispersive channels which 
may be fading is considered. In section 1 a VA receiver 
structure using whitened matched filtering is considered. In 
section 2 the Viterbi algorithm is described. Its perfor- 
mance is considered in section 3 and in section 4 a 
receiver structure for tropochannel applications using 
circulating filters is described. In section 5 some simula- 
tion results for tropo and telephone channels are given. 


4.1 VA Receiver Structure 


A receiving scheme 
decoding is shown in Fig 
matched filter, a symbol 


employing the concepts of Viterbi 
. 4.1. It consists of a whitened 
rate sampler and the Viterib&Hle coder. 
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For the case of finite inpulse response channel and 
white noise the output of a. filter h(-t) Hatched to the 
channel h(t) forms sufficient statistics for the estimation 
of the input sequence [Tan-68 J. Thus the matched filter 
outputs 


= J r(t) h(t-kl) dt 

_ co 

where r(t) is received signal, form sufficient statistics 
for the estimation of the input sequence a^. Since the 
received signal r(t) is 


r ( t ) = l a k , h(t-k'T) + n(t) 
k' 


Xu can be written as, after interchanging integration with 


summation 


00 CO 


X k = l a k J h(t-kT) h(t-k'T) dt + / li(t) h(t-kT) dt 


k 


, “k'i 


= l 


k 


a vi \_k’ + n k 


i 


where 

oo 

%_l > = / k(t-k'T) h(t-k'l) dt |k - 1:' f < n - 1 

•.oo 

_ q Otherwise, 

is the pulse auto -correlation coefficient of h(t) and there 
are 2n-l coefficients. Using the delay power series 
representations 
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x(D) = x Q + x^D + XgD 2 + ...... 

a(D) = a Q + a-jD + a^D 2 + ...... 

S(I» = S.( n _i) D” (n “ 1} + R. n+2 D- :1+2 + 
n 1 (D) i n;_ + rip + ny> 2 + ..... 


•• + B ^" 2 + Vl ^" 1 


¥ e can writ e 

x(D) = a(D) R(D) + n' (D) 

Here n' (D) is zero nean coloured Gaussian noise with, auto- 
correlation function a R(D) where a is spectral density 
of n(t). Due to the coloured noise, the filter output cannot 
he used in the Viterbi decoder and whitened ns.tch.ed filter 
which gives uncorrelated noise at its output is needed. 

Since R(D) = R(D~*^) t its roots occur in pairs of the 
form (jB, p"" 2 "). Collecting one root from each pair polynomials 
f(D) and f(D~ 2 ') \.ro formed such that 


R(D ) = f(D) f (D - ” 1 ) . 

Then we can write the coloured noise as 
n 1 (D) = n(D) f(D -1 ) 

since n* (D) then has variance a 2 f (D f (D) = cr 2 R(D) 
and zero mean Gaussian process is completely specified by its 


variance 
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Then , 

x(D) = a(D) f(D) f(D“ 1 ) +■ :i(D) f(D -1 ) 
or, 

z(D) = = a (D) f (D) + 11(D) . 

f (D ) 

in which, noise is white. 

The sequence z(D) can be obtained by sampling the 
outputs of the cascade, termed whitened MF, of a MF h(-t) 
and a transversal filter l/f(D“ 1 ). The whitened MF is 
realizable when f(D "**) has no roots on or inside the unit 
circle. 

That the sequence z(D) forms sufficient statistics for 
decisions about the sequence a(D) is clear from the fact 
that the sufficient statistics x(D) can be recovered by 
passing z(D) through the filter f(D - '*'). The Viterbi 
decoder then operates upon the sequence z(D) to generate 
Ml estimates of the sequence a(D). 

4 . 2 The Viterbi Algorit hm 

The input to the Viterbi decoder is 

z(D) = a(L) f(D) 4 - 11(D) 

= y(D) + 11(D) 

The decoder obtains ML estimate of the sequence a(D) by 
maximizing the log likelihood function 



56 


A' =' In p[z(D) / a ( D ) ] 
or equivalently 

A = In p[z(D) / y(D)] 

Since the noise samples are independent Gaussian distributed 
A = J In P|.z k /y k ] 

■where 

p[z 

¥e define 

k-1 

A k = I in Pl> k /y k ] 

3c- l = 0 11 

and 

6 A k = In p[z k /y k ] 

so that 

A k+1 = A k + 6 A k 

The log likelihood function is thus represented as a 
sun of independent increments. The following description 
of the algorithm illustrates how the maximization of the 
likelihood function is recursively achieved in Vitorbi 
decoding. 


k/y * ] = 8XPL 


T z k r_ y h . 
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The Viterhi algorithm can bo described either by 
finite state machine approach or by trellis diagram method. 

We follow the second approach where the lowpass equivalent 
filter h(t) of the channel is modelled as convolutional 
encoder . 

The filter h(t) can be visualized as convolutional 

encoder [G-ou-72] using its transversal filter representation 

shown in Pig. 4.2, The corresponding trellis for constraint 

length I = 3 is shown in Fig. 4.3. The trellis assumes a 

fixed periodic structure after a depth L. After this point, 

each of the 4 nodes (in general nodes) at time k are 

reached from two nodes corresponding to tine (k-l). Thus 

there are two paths corresponding to two dal a sequences 

entering a given node, j ; j = 1,2,3, .... 2 Ir " 1 . The 

Viterbi algorithm calculates the probabilities of transition 

to the node j, the incremental log likelihood functions 

6 A, .. , in the interval (k+1) from various nodes in tho 
^ 9 3 

interval k. These 6A V • are used in updating the log 
likelihood functions • of the path i j i = 1,2, entering 

9 J 

the node j. The algorithm eliminates all but the most likely 
path from further considerations by retaining the path with 
largest A^ Thus at time (k+1) there are 4 (in general 
2 1- ' 1 ') retained paths or survivors with log likelihood 


functions 
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Vitorbi [7it-71j has shown, that a 
than 5.8 times the constraint long 
results in only negligible degrade. 

4.5 x orfo cc_ Analysis 


decision delay of greater 
th is sufficient end 
tion of the performance. 


4.5.1 Probability of Error Performance 

'i'no Viterbi decoder is a non-linear decoder and its 
performance cannot be analyzed in closed form Certain 
bounds uii its performance have been obtained by Viterbi 
[Vit-'/l] and homey [for-72j. Forney has shown that at 
medium and high SNil' s the symbol error probability is over- 
bounded by 




who r o ^ is a snail constant, d n±n is the mininun energy 
of any n m-zero signal and is the spectral density of 
the noise and *(•) is the probability of error function. 

On the other hand , for any estimator the symbol error 
prob. is lower bounded by 


P(e) 


_ 




i 2o j 


where .<£ 
order of 


is another constant, typically within an 
magnitude of i^. Thus the probability of error 


is tightly bounded by 
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*‘L ; ^ ‘;iia /2a j - Jr '( £ ) i - vfhd. \1tJ 2o ] 

.ji i^ a - 2, • i>e ’iuito close, tii is result implies that with, 
vie. ibi . ... o j ... in,_, w.i ,.io3 g channels ISI docs not load to 
r.i.y : 1 ,, . ii i j- at uogradation in performance, and the- performance 
in ■.li.urt optimum. In fact, the structure is optimum for Ml 
cctinati -.i v f tiro entire transmitted sequence [For~72j. 

‘7*3*2 implementation of the Viterbi Decoder 

The* da cod or may be constructed in a completely parallel or 

oorivJL o. if L .urntioii or in a mixed configuration. In the 

pa r-.il el pr*. ..•cursin';; 2^ magnitude devices (adder + storage) 

, 1 

•arid 2“ comparators arc required. The fastest computation 
opcode arc I/i ilz. In serial configuration only one magnitude 
device and j.ie comparator arc needed but the fastest computa- 
ti.ii speeds are 2/1 Hz. In serial/ parallel configuration 
speeds n' 2 */ xT are required where x is a complexity factor 
which dop mci wi: tno proportion of serio.l to parallel 
circuitry . 

,-m,r data rates in the range of Mbits, sec. the hardware 
ruaiircmcnts of a Vitorbi decoder aro unrealistic with 


pr.;Ou 

lit £ 

cy 

t ochn-1 

inv ..1 

vet. 


The data 

t olep 

ilwUC 

- 

iiannols i 

on ou 

! i c, 

i.a: 

miclc . 

tioli 

aris 

cs 

duo to 


these voiceband channels limita- 
ISI. For these channels 1 may 
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be aw 1 ' 0 • r - G 4u or 50. Though high speeds of operation 
arc not required a largo number of magnitude devices and 
co..-.?ar:\tvns arc needed and the system is not econonical. 

4.4 y\. A-ccivor Using Circulating Filters 

An mi tinned earlier tne input to the Viterbi decoder is 
obtained by sampling the output of a whitened Hatched filter 
w(t) whose input is the received signal r(t). 


r = / r(t) w(t~k!) dt. 

_ 00 

The* white.. od matched filter is show/: to bo a cascade of a 
filter . sat died li(-t) to the* channel h( t) and a transversal 
filter, whose tap gains are given by the sequence 1/f (if 1 ) 
where 1(D ) is defined as 


Kir 1 ) = t + qc" 1 + q d 2 + 


+ ~n-l 1 


-(n-1) 


it is p minted out in j_For— 72j that cno cransversal filter 
is realisable when f (D has no roots on or inside unit 
circle. knowledge about channel inpulse response is needed 
in both the Hatched filter and the transversal filter. 

In Fig. 4,4 a VjL receiver scheme which uses circula- 
ting filter in estimating- the channel impulse response xs 
shown. The channel estimates obtained by the circulating 
filters arc used in whitened matched filter in setting 
the tap gains of the transversal filter, and also in matched 
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.. iltjring die received, signal, a, decision delay t of 
6 to 7 tines the constraint length, is involved in the 
receiver decisions the channel estimates arc available after 
a a clay ■> f x sec. and the assumption to be nad© is that tile 
cnar oh-riotics do not vary significantly during this period, 
fhe as sumption is valid for telephone charnels as. well as 
for high data rate transmission over troposcatter channels. 
The waxte-ned matched filter output signal is then operated 
upon by the Vltorbi decoder to obtain the ML estimates of 
the transmit tod symbols, 

/or slowly fading channels it nay not be required to 
adjust the tap gains of the transversal filter eveiy I sec. 
(adjusting tap-gains involves complicated complex calcula- 
tions), The tap gains can be adjusted every pi bits where p 
nay bo a few hundred depending on data rate and fading 
characteristics of the channel. This reduces the speed 
requirements on hardware needed. In simulation it reduces 
Compulation Lime. 

4 , 5 oi. iulati on. Res ults 

Details of simulation strategy followed are given in 
the Appendix 1. Due to largo amount of computation time 
involved only one 'run' of the program could be obtained in 
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which no errors were encountered for 6:xl0^ data bits. The 
Slid was 6 dB and diversity reception was not used. The 
channel was an unknown random dispersive channel. 

Some results were also obtained for telephone channel 
applications. These results along with the samples of 
channel response used are given in Table 5. 

Table 5 

Performance of Viterbi decoding in telephone 
channel equalization 

g~o.__of errors in 10^ bits^ 

26 
2 

1 


ShR(dB) 

5 

10 

15 


* Results averaged over 5 different noise sequences for 
same input data string 

Telephone channel characteristic - 


40.1, -8. 3, “16. 6, 3.0, 19.6, “0.5 



CHAPTER 5 


CONCLUSIONS 

In this thesis we developed an adaptive decision- 
directed channel estimation technique using circulating 
filters and studied its p erf o rum ce in two different 
equalization schemes for troposcatter channels. The equali- 
zation schemes selected, the decision feedback equalization 
and the Viterbi decoding, are representative examples of two 
different classes of nonlinear equalization techniques. 

The first technique is the best among the techniques practical 
with present day technology and the second one is an optimum 
technique which is inplementable too. 

The channel estimators using circulating filters have 
the desirable property of suppressing noise. The theoretical 
derivation indicates that the circulating filters can 
suppress the noise in channel estimates by more than 12 dB. 
This substantial noise suppression is verified by simulation 
results. The capabilities of the scheme are further mani- 
fested by the simulation results that small time jitters and 
error propagation due to small bursts of errors (upto 5 
errors) have no apparent effects on channel estimates. 

A DEE receiver scheme for troposcatter channels which 
has no forward filter is developed by taking advantage of 



68 


highly skewed nature of the delay power spectrum of the 
channel. The receiver thus does not have any performance . 
loss due to inadequately synthesized forward filter [c.f.. 
Monsen's Structure, Mon-74]. Simulation results for such a 
DfE receiver using circulating- filters for channel estimation, 
though inadequate for making any definite observations, suggest 
that the system has a quite low probability of error. The 
results without diversity compare well with the theoretical 
performance bounds. Simulated error propagation effects are 
found to cause negligible performance degradation for small 
bursts of errors, furthermore, small time jitters of + T/4 
are seen to have insignif icant effect on error rate. 

The simulation results for Yiterbi decoding using circu- 
lating filters are quite insufficient and no concrete 
conclusions about the receiver performance for t'roposcatter 
channels can be derived escept that, possibly, the decoder is 
capable of performing well in an unknown channel condition. 

Unfortunately, due to large amount of computation times 
involved sufficient data for making conclusive observations 
about the system performance could not be obtained. Certain 
assumptions had to be made to reduce the computation times 
required. 'Thus, in simulation of D .• E receiver error propa- 
gation effects were ignored in the computation of the average 
probability of error. That this is a valid assumption for low 
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error-rate systems is justified by simulation results on 
error propagation effects. In Yiterbi decoding it was 
assumed that the adaption of the whitened matched filter once 
in a thousand bits is sufficient, for high data rates of 
transmission over slowly fading channels lime troposcatter 
channels this is a reasonable assumption. 

The simulation programs not requiring much computation 
time were run on IBM-7044 compucer at 1. 1. 1, iianpur,, 'The 
simulation results on average probability of error were 
obtained on DB0-10 system at TIPR* Bombay. It took 9 hours of 
computation time on 1)1)0-10 system to obtain results for 
DYE receiver. Another one hour was taken by VA receiver pro- 
gram. Total cost incurred in simulation on DBO-lp system 
was six thousand rupees. 

One clear possibility of further work on the problem 
considered in this thesis is to obtain simulation fyata for 
more channel realizations and larger data sequences? and to 
study the adaptive behaviour of the receivers by varying the 
delay power spectrum profile. Since time synchronisation is 
extremely difficult in tropo channel, applications it is 
desirable to study in detail the effects of time- jitter 
on the performance of the equalizers. 
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un the theoretical side work can he done on developing 
a DJ?3 scheme which utilizes all the received signal energy 
without suffering any performance degradation due to 
inefficient forward filtering. 
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APP JED IX 1 


PERFORMANCE S MUTATION 

Phe objective of the computer simulation was to study 
the behaviour of an adaptive channel estimator using cir- 
culating filters and to evaluate for troposcatter channels 
the performance of a decision feedback: equalizer and a 
Viterbi decoder with circulating filters for channel estima- 
tion. Sampled versions of continuous-time systems were 
considered. A trade-off between number of samples per baud 
and computation time required was found necessary. The 
simulation was carried out with 4 (M) samples/baud. 

As described in Chapter 2, a troposcatter channel is 
completely described by specification of its delay power 
spectrum. The channel impulse response is selected from an 
ensemble defined by an approximation to the delay power spe- 
ctrumj Q(t) = a xe~ . The complex zero-mean uncorrelated 
Gaussian samples with variance equal to the delay power 
spectrum at the sampling instants constitute the channel 
impulse response. The channel impulse response is finite 
with a length of 3(N) bauds. 

The data input to the channel is a pseudo-random binary 
sequence with symbol separation T normalized to unity. The 
sequence is obtained by generating a random number using the 



subprogram RUDY I (RAF on DEG-10 system). If the random 
number is greater than 0.5 the input symbol is taken as 
+1, otherwise -1. 

The channel output is the convolution of the transmitted 

data sequence and the channel pulse response. The channel 

pulse response is obtained by convolving channel impulse 

response with a unit pulse of normalized baud duration. The 

received signal is then obtained by adding a zeromean white 

2 

Gaussian noise sample with variance a to the channel output. 
The Gaussian noise sample n^ is generated by using the Box- 
Muller formula 

n^ = o[-21n (R- l )] t cos(2tcR 2 ) 

where R^ and R 2 are uncorrelated uniformly distributed 
random numbers in the range 0 and 1. The noise variance is 
adjusted to fix the SM of the received signal. 

The steps involved in obtaining the received signal are 
given in flow-chart 1. 

In the following section, the estimator behaviour in 
extracting the channel information from the received signal 
is studied. In the subsequent sections the decision 
feedback equalizer and the Yiterbi decoding schemes with 
circulating filters are simulated and their performance 
evaluated. 



A 1.1 Channel Estimation Scheme 


The estimator is used to determine the channel pulse 
response. For a channel impulse response of duration 
^(=3), two sets of (N+l) circulating filters,* one for the 
quadrature component and the other for in phase component of 
the received signal, are needed. 

The simulation steps are shown in Flow chart 2. 

The sample at time T/4 of a given chip, h°(t) is estima- 
ted in the jth circulating filter in the first iteration 
and the sample at time T/2 in the second iteration and so on. 
In estimating the samples of the chip h?(t) the interference 
contribution due to other chips is subtracted from the 
received signal sample using past decisions and channel 
estimates in the preceding iteration. Correlation with the 

symbol b, . gives a noisy channel estimate . The circulating 
k — 3 

filter averages these estimates in successive iterations by 
adding its input signal to the filter output signal in the 
preceding iteration weighted by feedback gain f . . Th^ 

J 

channel estimates are given by circulating filter output 
weighted by (1-f.). 

<J 

The noise suppression due to the circulating filter 
was studied by computing the difference in SHR's at the 
filter output and input. Effects of decision error on 
channel estimates were investigated by deliberately 



introducing decision errors. The time-jitter effects were 
evaluated by shifting the received signal in steps of + T/4. 

A 1.2 Decision Feedback Equalization 

The simulation steps are shown in Flow chart 3. The 
samples of the received signal are generated as described 
earlier. The channel estimates are obtained from the 
received signal using receiver decisions. The channel 
response estimates and past decisions of the receiver are 
used in computing the past symbol interference. The signal 
samples are then multiplied by the corresponding samples of 
the channel response in the first T sec. The resulting signals 
are summed up to form sampled output of the matched filter.. 

If it is greater than 0 the decoded bit is -i- 1 otherwise -1. 

The receiver decision is checked against the transmitted bit 
to detect decision errors. 

The steps are repeated over a large number of input 
data bits and also for different realizations till sufficient 
data for computation of the probability of error is obtained. 

The performance was studied mainly in the context or 
evaluating the average probability of error and studying 
the effects of .error propagation and time jitter. The last 
two effects were studied by deliberately introducing decision 
errors and time-jitters in the received signal respectively 



The computation of error probability required a very large 
amount of the computation time. To reduce this time certain 
assumptions were made. Thus, it was assumed that the receiver 
decisions are essentially correct, and the error propagation 
was not considered. Another assumption was that in diversity 
reception there is a very low probability of decision error 
when a correct decision is obtained in one branch. This 
implies that the need for considering all diversity branches 
arises only when an error occurs on a single branch. 

A 1.3 Viterbi Decoding 

The flow diagram of simulation is shown in Plow chart 4. 
The channel impulse response estimates are used in the 
computation of the whitening sequence f(D). First, the 
sequence R(D), autocorrelation function of impulse response 
h(t), is obtained by convolution of h(t) with itself. The 
roots of the polynomial equation R(D) are then computed and 
the sequence f(D) is formed by collecting those roots which 
lie outside the unit circle,- i.e. which have magnitude greater 
than unity. 

Since computation of the sequence f(D) involves compli- 
cated complex calculations, the sequence is updated only 
once in a thousand bits. It is assumed that the data 
transmission is at sufficiently high speed over a slowly 
fading channel, so that the channel characteristics do not 



change significantly in the thousand bits duration. 

The input to the Viterbi decoder is then obtained by 
convolving the input data sequence with the sequence f(D) 
and adding a white Gaussian noise sample to it. 

As described in Chapter 4 , the Viterbi algorithm pro- 
ceeds recursively in two steps. It calculates the log like- 
lihood functions for all allowable paths leading to a given 
node in kth iteration by calculating the associated probabi 
lity of transition and adding it to the log likelihood of tb 
path in (k-l)" fc * 1 interval. In the second step it chooses the 
path with maximum log likelihood as the survivor. As deri- 
ved in [Vit-71], for BPSK transmission, the probability of 
transition is given by the product of the code symbol and 
received sequence. The path histories, which are in the 
form of sequences of +1 and -1 are stored in their decimal 
equivalent form. The decoding delay is kept at 7 times the 
channel impulse response spread and the bit corresponding 
to the first node is taken as the receiver delay. 
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